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Abstract 

A  theoretical  investigation  is  undertaken  into  the  creep 
instability  of  two  common  shells  of  revolution,  the  infinitely 
long  cylinder  and  the  complete  sphere-   The  linear  and  non- 
linear visco-elastic  material  idealizations  are  used  in  the 
various  theoretical  procedures.   The  critical  buckling  time 
is  greatly  influenced  by  the  parameter  a/h  ratio,  magnitude 
of  initial  imperfections  and  material  selection.   The  spherical 
shell  has  a  lifetime  longer  than  the  cylindrical  shell  having 
the  same  a/h  ratio,  initial  imperfections  and  material  con- 
stants, while  the  critical  mode  numbers  are  the  same. 
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NOTATION 


w  transverse  displacement,  positive  outwards. 

u  displacement  in  n   direction. 

1 

v  displacement  in  r\      direction. 

2 

M  ,M  ,M   ,M  moments  per  unit  length  as  defined  in  Figure  2. 
"   2   21   12 

,N  ,N   ,N  forces  per  unit  length  as  defined  in  Figure  2. 
2   21   12 

,e   ,e   ,e  total  strains. 
1   22   21   12 

,P  ,P  external  loads  per  unit  area  as  defined  in  Figure  2, 
2   3 

,x  ,x  orthogonal  coordinates  as  defined  in  Figure  1. 
2   3 

,Q  transverse  shear  forces  per  unit  length  as  de- 


ll 


2 


fined  in  Figure  2. 


,R  Principle  radii  of  curvature  of  a  shell  and 

radii  of  curvature  of  the  coordinate  curves 


n  ,n  . 
1   2 

a  cylindrical  shell  mean  radius. 

a*  out-of -roundness  parameter  defined  by  equa- 

tion (4.19). 

h  cylindrical  shell  mean  thickness. 

S  arc  length. 

A,B  defined  in  equations  (2.5  b-c) . 

2  distance  measured  along  a  normal  from  the 

middle  surface  of  a  shell. 

E  Youngs  Modulus  of  Elasticity. 

C,K  material  viscous  coefficients. 

M  creep  exponent  defined  by  equation  (3.4). 

R  spherical  shell  mean  radius 

R*  out-of- roundness  parameter  of  spherical  shell, 
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NOTATION  (cont.) 


H  spherical  shell  mean  thickness. 

A  amplification  function  (cylinder) . 

n 

t  time . 

t  critical  time, 

cr 

B  amplification  function  (sphere) . 

mn 

C  /D  maximum  initial  imperfection;  cylinder/sphere, 

n   mn  c 

n  harmonic  in  circumferential  (6)  direction. 

m  harmonic  in  circumpolar  (cf>)  direction. 

H  ,n  lines  of  principle  curvature. 

1   2 

k  ,K  change  of  curvature. 

1   2 

k  half  total  change  of  curvature  in  twist. 

12 
*   * 

K        K 

I      2  principle  curvatures 


R 

1  2 

i 

£  ,£  membrane  strains. 

1   2 

Y  shear  strain. 

0 

u  poisson's  ratio. 

o*  effective  stress  defined  in  equation  (3.5b). 

4   ,     4 
a  —  Ch  or  —  CH. 

6,x  circumferential  and  axial  coordinates  in 

cylindrical  shell. 

4>,6  circumpolar  and  circumferential  coordinates  in 

spherical  shell. 

4        4 
ft  —  kh  or  —  kH. 

Ch3    ai3 

«  -?"  or  T"  ' 

Kh3     KH3 
C  -7—  or  — — -  . 
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NOTATION  (cont.) 


X  -aaw/C  (time  constant  linear  viscous  cylinder) . 

X  -afiw/C  (time  constant  non-linear  viscous 

cylinder) . 

$  -3Raw/2£  (time  constant  linear  viscous  sphere) . 

A  -3RTiw/2C  (time  constant  non-linear  viscous 

sphere) . 

t  Xt  (non-dimensional  time  linear  viscous  cylinder) 

t  Xt  (non-dimensional  time  non-linear  viscous 

Q 

cylinder) . 

t  $t  (non-dimensional  time  linear  viscous  sphere) . 

s  ,      c 

t  At  (non-dimensional  time  non-linear  viscous 

sphere) . 

A  n(n  +  1) 

n 

V  defined  by  equation    (5.12c). 

ij>  o*/a* 

S.  . 

3-D 


0 

stress  deviator. 

initial  value  of  ( 

). 

3(   ) 

9i 

3(   ) 

3t 
dominant  part  of  (   )  . 
perturbation  part  of  (   ) 
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1.   INTRODUCTION 

1. 1  Preface 

The  problems  of  creep  deformation  and  buckling  of  structural 
members  while  often  addressed  in  aerospace  and  chemical  engineering 
have,  until  recently,  been  aphonic  design  considerations  in  ocean 
engineering.   The  rationale  being  that  the  ocean,  itself,  offers  an 
environment  very  unaccommodating  to  creep  with  its  relatively  low  tem- 
peratures.  Further,  the  general  use  of  creep  resistive  materials  for 
relatively  short  periods  of  loading  made  creep  insignificant  in  the  de- 
sign of  ocean  shell  structures.   Recently,  the  escalating  demand  man 
and  his  society  have  placed  on  the  ocean  in  pursuit  of  its  understand- 
ing and  resources  has  resulted  in  a  rapid  advancement  in  ocean  technology. 
This  new  technology  has  generated  an  increased  interest  in  the  creep  of 
shell  structures.   This  manifestation  of  interest  has  been  stimulated 
by  many  events,  the  greater  use  of  creep  prone  materials,  aluminum  and 
titanium,  the  longer  duration  structures  are  required  to  withstand  com- 
pelling hydrostatic  pressures,  and  the  need  to  predict  volumetric  changes 
due  to  creep  of  shells  which  have  stringent  volumetric  design  constraints, 
such  as  the  neutrally  buoyant  and  compressible,  free  floating,  current 
measuring  SOFAR  floats  of  the  Woods  Hole  Oceanographic  Institution. 

In  the  past,  the  ocean  engineer  was  perplexed  with  a  two  part 
dilemma  when  asked  to  predict  the  creep  performance  of  his  design,  first, 
find  an  acceptable  analytical  model  for  creep  deformation  of  shells  and 
second,  locate  accurate  material  constants  for  low  temperature  creep. 
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The  purpose   of   this  paper  is   to  solve   the    first   dilemma  of 
developing   a  simple   analytic  technique   to  predict   creep   deformation   of 
two   common   shells   of   revolution,    the    cylinder   and  complete   sphere,    suit- 
able   to  oceanic   applications    and  which  will  be   useful   in  preliminary 
stages   of  design   and  may  be   adequate    for   final   design.      Finally,    to   dis- 
cover if   creep  buckling  is   a  phenomenon   that  needs   to  be    actively   con- 
sidered in   the   design   of   long  term  ocean   shell   structures. 
1.2      Historical   Review 

A  literary   search   of   the   existing   literature   reveals   that   few 
published  works   have  been  transcribed  which   are   directly   applicable   to 
conditions   experienced  in   the   ocean.      Most  papers   deal  with   columns, 
plates   or   cylindrical   shells    loaded  in   axial   compression.      A  number  of 
investigators  have    attempted  solutions    for  the    creep  behavior  of  thin 
cylindrical   shells    loaded  under  uniform  external  pressure.      To   the   best 
knowledge   of  this    author,    no  paper  has   been  published  dealing  with   the 
creep  performance   of   a   complete   spherical   shell   loaded   under  external 
pressure. 

The   author's   dealing  with    cylindrical  shells   use   a  variety 
of  simplifications,   but   do  have    an   assemblage   of   assumptions    common   to 
all   their   analyses.      These   similar   assumptions    include: 

(a)  thin   shell   assumption,    so   that   the    radial   stress    compo- 
nent  is   negligible   in   comparison   to  the   in-plane   stress    components 

(b)  only   the   effect   of   secondary   creep,    according   to   Odquist's 
law    (1937)    are   considered 

(c)  buckling  is   not   defined  in   the  bifurcation   sense,    but 
rather   as    a  state  of  progressive   accelerated  growth   rate   of  the   displacement 
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field,    creating  unfavorable   strain   distribution   and   leading   to   full 
structural    failure. 

(d)  initial   imperfections    of   a  shell   cause   eccentricity   of 
loading,    thus    serving   as    a   lever   arm   for   compressive   in-plane    forces, 
which   in   turn   increases    the  bending  moment,   which   in   turn,    increases   the 
rate   of   curvature    change 

(e)  an   assumed  solution   of   an  infinite   trigonometric   series, 
either   a   complete   or  partial   Fourier   series. 

Hoff,    et    al    (1959)    used  a  sandwich  model,    in  which   the   outer 
skins   support  membrane    and  bending   stresses    and  the    core   supports   shear 
without   deformation.      The   sandwich   model   avoids    complications    in   the 
superposition   of   compressive    and  bending  stresses   due   to  the   non- 
linearity   of   the    creep   law.      Further,    a  quasi-elliptical  initial   imper- 
fection was    assumed.      Samuelson    (1970)    developed   a  numerical   solution 
for  non-uniform  external    loads.      Sankaranarayan    (1969)    generated  a  solu- 
tion   for   creep   deformation   of   cylindrical   shells    under   combined   lateral 
and   axial  pressures    for   a  material   that   adheres    to   the   Tresca  criterion. 
Berman,    et   al    (1974)    presented   a  numerical   analysis    for   a   cylindrical 
shell   infinitely    long,    so  that,    a  state   of   generalized  plane    strain   exists. 
As   did  Hoff,    Berman   assumed  that   the   initial   imperfection   is   quasi-elliptical. 

Penny  and  Marriott  (1969)  have  carried  out  the  only  known,  but 
rather  special  analysis  of  creep  deformation  of  a  spherical  shell.  They 
explored  snap   through   of   a  spherical  boss  point    loaded. 

All   analyses   emphasize   the    fact   that   initial   imperfections 
play   a  most   important   role   in   creep  buckling.      That,    buckling  stems 
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from  the   accelerated  growth   of   these   initial   imperfections   in    an   otherwise 
uniform  initial   displacement    field. 
1. 3      Method  of  Analysis 

The   conviction   that   creep  buckling  stems    from  the   growth   of 
small   imperfections    in   the   otherwise   uniform  initial   displacement   field 
suggests   the   use   of   a  perturbation  method  of  solution   as    first   developed 
by  Abrahamson   and  Goodier    (1952)    for  dynamic  plastic  problems.      The 
method   assumes   an   axisymmetric   dominant  solution   in  which  buckling   does 
not   occur   and   a  perturbation   solution   in  which  the   influence   of  bending 
is    significantly   greater   than   the   in-plane    forces   thereby   giving   rise 
to   inextensible  buckling. 

Due   to  the    axisymmetry   of   the   dominant   displacement    field  of 
the   shell  of   revolution,    only   the   transverse   displacement   field  consists 
of  a   dominant   and  perturbation  part,    while   the   in-plane   displacements 
have   only   a  perturbation  part   so  that 


w(n    ,    n    ,    t)    =  w(t)    +  w'  (n    ,   n    ,    t)  (1.1a) 

12  12 


u(n  ,  n  ,  t)  =  u'(n  ,  n  ,  t)  d.ib) 

12  12 


and 


v(n    ,   n    ,    t)    =  v'  (n   ,   r\   ,   t)  (1.1c) 

12  12 

respectively,  with  perturbation  parts  delineated  by  primes  and  dominant 
parts  by  bars. 

In  turn,  to  remain  consistent,  the  shear  strain  and  twisting 
curvature  change  have  only  perturbation  parts  while  the  membrane  strains 
and  bending  curvature  changes  have  both  dominant  and  perturbation  parts. 
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Similarly,    the   shear   forces   and  twisting  moments   have    just  perturba- 
tion parts  while  normal    forces    and  bending  moments   have   both   dominant 
and  perturbation  parts 


N    (Ti    f  T\    ,  t)    =  N    (t)    +  N'     (n    ,    n    ,    t)  (1.2a) 

112  1  112 

n  (n  ,  n  ,  t)  =  n  (t)  +  n'    (n  ,  n  ,  t)  (i.2b) 

2       12  2  2       12 

n     (n  ,  n  ,  t)  =  n«      (n   ,  ti  ,  t)   =  n     (n  ,  n   ,  t)         (i.2c) 

1212  1212  2112 

m  (n   ,  n   ,  t)   =  m  (t)   +  w    (n  ,  n   ,  t)  (i.2d) 

X      1         2  1  112 

m  (n   ,  n   ,  t)   =  m  (t)   +  m"    (n   ,  n   ,  t)  -       (i.2e) 

2       12  2  2       12 


and 


m     (n  ,  n  ,  t)  =  m'     (n  ,  n  ,  t)  =  -m     (n  ,  n  ,  t)       (i.2f) 

12       1         2  12       1         2  21       1         2 


The  material   comprising   the   shell   is   isotropic,    incompressible 
and  homogeneous.      The   initial   elastic   response    is    considered  to  be    infini- 
tesimal.     No   initial  shape   of  imperfections    is    inferred.      Further,    the 
assemblage   of   assumptions    common   to  past   analysis,    (a)    -    (e)    on  page  11 
is    assumed. 

The    analysis   will  present  non-linear  equilibrium,    strain- 
curvature   relations    and   constitutive   equations    for  linear  and  non- 
linear visco-elastic  materials    for   an   arbitrary   shaped   shell.      The 
reduction   of   these   relations   will  be    carried  out    for  the   infinite    cylinder 
and  the    complete    sphere. 
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2.   BASIC  RELATIONS  FOR  THIN  SHELLS 

2. 1  Introduction 

The  succeeding  relations  were  developed  by  V7alters  and  Jones 
(1972) .   They  were  derived  with  the  aid  of  differential  geometry  as 
discussed  in  Kraus  (1967)  and  are  valid  for  any  thin  arbitrary  shell 
which  has  lines  of  principle  curvature  aligned  with  the  parametric 
lines. 

2 .2  Shell  Coordinates  -  Differential  Geometry 

The  coordinates  of  point  P  on  the  middle  surface  of  a  shell, 

the  reference  surface,  are  defined  in  Figure  1.   r  is  the  position 

vector  which  is  a  function  of  parametric  lines  n  ,  n  .   n  is  the  unit 

1    2 

normal  vector.   The  position  vector  can  be  expressed  as 

r=xi   +   xj   +   xk  (2.1) 

1  2  3 

/S   /\  /\ 

where   i,    j,    and  k   are   the   unit  vectors   in   the   x    ,    x    ,    x     directions, 

1         2         3 

respectively. 

As  demonstrated  in  Kraus  (1967) ,  the  magnitude  of  the  prin- 
ciple curvatures  are  the  roots  of  the  following  quadratic  equation 

H2K2  -  (EN-2MF-GL)K  +  (LN-M2)  =0        (2.2) 
where 

E  =  r ,  -r,  (2.3a) 

1    1 

F  =  r,   r,  (2.3b) 

1    2 

(equal  to  zero  when  n   is  orthogonal  to  n  ) 

1  2 
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G  =  r,  «r,  (2.3c) 

2    2 

H  =  (EG-F2)^  (2.3d) 

S\  — 

L  =  n*r,  (2.3e) 

11 

/\  — 
M  =  n-r,  (2.3f) 

12 

(equal  to  zero  when  r|   and  n   coincide  with  the  lines  of  curvature) 

1      2 

N  =  "n^r,  (2.3g) 

22 

with  —   — 

r,  Xr, 

/\     *    12 
n  = >    (2.4) 

Equations  5(a-d)  are  referred  to  as  the  first  fundamental  magnitudes 

while  equations  5(e-g)  are  called  the  second  fundamental  magnitudes. 

If  n   and  n   are  orthogonal,  then  F  equals  zero,  the 
1       2 

elemental  area  can  be  expressed  as 

ds2  =  A2d2n   +  B2d2n  (2.5a) 


where 


and 


A  =  (E)  2  (2.5b) 


B  =  (G)"5  (2.5c) 


2. 3   Strain-Curvature  Relations 

\     A  tensor  formulation  and  the  usual  assumptions  of  shell 
theory,  (eg.  Kraus  (1967))  have  been  employed  by  Walters  and  Jones  (1972) 
to  show  that  the  strains  in  the  thin  shell  may  be  written 
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and 


e        =   e      +  3k  (2.6a) 

11  1  1 


e        =   e      +   2k  (2.6b) 

22  2  2 


e        +   e        =   y      +   2Zk  (2.6c) 

12  21  0  12 


where   the   membrane   strains    are 


u,  A,    v  (w,    )2 

1  2  w  1 

£      -  +  — —  +  —  +  (2.7a) 

!  A  AB  R  2ft2 

v,  B,    u  (w,    )2 

2  1  2 

w 

e      =  — — -  +    ■-_     +  — -  + >        (2.7b) 

2  B  AB  ^  2r2 


and 

w,   w, 

A,u,  B,v, 


\  =  b<a',   +  a'i*,   +-isr-  (2-7o> 

2  1 


and  the   curvature   changes    are 


w,  A,  w. 


1    u  1  2     v  *■ 

k      =  -(— )  +  ( )  (2.7d) 

I  A    R  A       ,  AB     R  B 

1  1  2 

w,  B ,  w, 

l     v  2  >      u  ' 

K2  -  f(it  -  — >  +  Tr'ir  -  — »  (2-7e) 

2  1 
and 

^■^"V.*^  ■"■?».       (2-7f> 


Note   that   y      is   the   total   shear  strain  while   k        is   one-half  the   total 
0  12 

twist   curvature    change 

The    Love -Kir chhoff   assumption  has   been   invoked   in   the 
derivation   of  these   equations. 
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2 . 4      Equilibrium  Equations 

The   typical   shell   differential  element  illustrated  in   Figure   2 
will  have    a  rate   of  internal   energy   dissipation. 


W.    ^   =    /In  e 
mt       J      i    i 

sL  J 


+  Ne      +  N      y      +MK 
2    2  12    0  11 


+   M  K      +   2M 
2    2 


K       [ds 
12    12J 


(2.8) 


where   s   is    the   total   un deformed  area  of  the   middle   surface   of   the 
arbitrary   thin   shell.      Now,    equating   the    rate   of   internal  energy  dis- 
sipation  to   the   rate  of  external  work    (i.e.,    using  the  principle   of 
virtual   velocities) ,    a   consistent   set   of  equations   of   equilibrium  is 
obtained. 


ABp      -    N    B,       +     (N      A)  ,       +    (N    B)  ,       +    N      A, 
1  2       1  12         2  11  12       2 


(MB),  M   B,  (MA),  M      A, 

11  2       1  12  2  12       2 


R 


R 


R 


R 


=   0  (2.9a) 


and 


ABp      -    N    A,       +     (N      B) ,       +     (N    A) ,       +   N      B, 
2  12  12         1  2         2  12       1 


(MA),            M   A,             (MB),             MB, 
2          2             12               12          1             12       1 
+ +   +   


R 


R 


R 


N    AB         N    AB         /N    Bw, 
1  2/1 

ABp + 

3  R  R 

1  2 


A 


+     (N      w,     )  ,       +     (N      w,     )  ,       + 
12       2       1  12       1       2 


R 


rN   Aw, 
2         2 


(M    B)  , 
1  1 

A 


=   0  (2.9b) 


■18— 


+  2M 

12,  12 

Inertia  terms  are  disregarded  since  creep  is  assumed  to  be  a  very  pro- 
tracted deformation  process 

3.   CONSTITUTIVE  RELATIONS 


3. 1   Linear  Visco-elastic  Material 

The  appropriate  model  to  simulate  secondary  creep  during  the 
application  of  a  step  function  load  is  the  Maxwell  model  illustrated 
in  Figure  3a.   The  strain  at  all  times  is  the  sum  of  the  elastic  strain 
component,  and  viscous  strain  component.   The  stress  in  the  viscous 
and  elastic  elements  are  always  equal.   Thus,  from  the  equation  of 
equilibrium  the  following  expression  for  strain  rate  is  found, 
(eg.,  Drucker  1967) 

;-i  +  § 

During  application  of  the  step  function  load,  the  motion  of 
the  dashpot  is  negligible.   Hence,  the  initial  strain  is  the  elastic 
stretch  of  the  spring.   As  time  procedes,  the  model  will  continue  to 
creep  steadily  at  the  rate  of  o/C  as  shown  in  Figure  3b. 

The  relation  for  the  linear  Maxwell  model  can  be  written 
in  tensor  notation  as,  (eg.,  Drucker  1967) 

ij    E  ij     E     3  ij   2  C 
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s 

• 

1+V-              3      ij 
=  s.  .    +  — J- 

e .  . 

i: 

E      ij         2      C 

where 

a   =  a   +  a   +  a 
Kl<     11     22     33 

5.  .  =  a.  . —6.  . 

6 .  .  =  kronecker  delta 

and  E,  C  are  material  elastic  and  viscous  constants 

respectively.   Since  superposition  of  hydrostatic  stresses  is  assumed 
to  cause  no  further  straining, equation  (3.2a)  reduces  to 

(3.2b) 

Inverting  equation  (3.2b)  and  performing  the  proper  intergration  reveals 
the  linear  viscous  constitutive  equations,  after  the  initial  elastic 
response,  that  is,  disregarding  the  elastic  behavior,  for  a  thin  shell  as 

N   =  a(e   +  •J-  e  )  (3.3a) 

1  1    2   2 

N  =  a(e  +  ^  e    )  (3.3b) 

2  2    2   1 

N    =   N    =  -r  Y  (3.3c) 

12     21    3   0 

M   =  £(K   +  i  k  )  (3.3d) 

1       1    2   2 

M  =  £(k   +  -  k  )  (3.3e) 


and 


where 


and 


2       2    2   1 


M   =  -M   =  £<  (3.3f) 

12      21      12 


a  =  |ch 


Ch3 


-20- 


3. 2   Non-Linear  Visco-Elastic  Material 

Secondary  creep  for  the  non-linear  visco-elastic  material 
can  be,  as  in  the  linear  case,  represented  by  a  Maxwell  model  where 

e--+  (E)  (3.4) 

This   model   does   not   simulate  primary   creep   and  has    a  poor 
response   to  variable   stress    conditions. 

Assuming  that   the    creep   strains    are   dependent   on   loading 
history,    the  Prandtl-Reuss    Incremental   relations    can  be   used   for  the 
non-elastic  part 

3   de* 
de.  .   =  ir^rS.  .  (3.5a) 

i]        2     a*    xj 

where 

a*  and  e*  are  effective  stress  and  strain  based  on  Von 
Mises  criterion 
where 

a*  =  -  [(a  -a  ) 2  +  (a  -a   )2  +  (a  -a   )2 


\(a     -a 
L  11 


/2     L  11   22        22   33        33   11 

+  6(t   2+t   2+t   2) I  ^   (3.5b) 
12    23    31  J 


e*  =  —  |~(e   -e   )2  +  (e   -e   )2  +  (e   -e   )2 
/J  L   1  1   22        22   33        33   11 

+  6(y   2+T   2+Y   2)  r   (3.5c) 
12    23    31   J 

Norton's  Law  for  secondary  creep  under  a  uniform  constant  load 

is  expressed  as 

*M 
de*  =  ?—-  dt,  (3.6) 
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Hence,  one  can  substitute  (3.6)  into  equation  (3.5)  and  find 

,   *M-1 

e..  =  IS- s..  (3.7) 

ij    2   K     ij 

which  is  Odquist's  Law  for  secondary  creep,  (eg.,  Odquist  (1956)). 

By  assuming  that  the  effective  stress  remains  constant  across 
the  shell  thickness,  equation  (3.7)  can  be  inverted  and  intergrated  to 
yield  the  non-linear  constitutive  relations  while  neglecting  the 
elastic  effects  (Penney  and  Mariott  1971) 

*(1-M)  •      1  • 
N  =  na  v   w  (e   +  ±   e  )  (3.8a) 

1  11    z      22 

*(1-M)   •       1  • 

N  =  Ua    v    ; (e   +  ±  e   )  (3.8b) 

2  22    2   11 

o*(l-M) 

N    =  N    = y  (3.8c) 

12     21       3       0 

*(1-M)  •     1  • 

M  -  r,a    K        '  (k  +  =-  k     )  (3.8d) 

1  11   2   22 

*(1-M)  •     1  • 
M  =  r,a    v    ;  (k  +  ~   k   )  (3.8e) 

2  22   2   11 

* (1-M) • 
M   =  -m   =  Ca      k  (3.8f) 


where 


12      21  12 


4 

n  =  —  Kh 


Kh3 


Note,  that  if  M  equals  one,  the  linear  case,  these  constitutive 
relations  reduce  to  the  linear  relations,  equation  (3.3). 
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4.   INFINITELY  LONG  CYLINDRICAL  SHELL 

4. 1   Introduction 

The  objective  of  the  following  section  is  to  reduce  the 
general  equations  for  an  arbitrary  shaped  shell  to  the  particular  case 
of  an  infinitely  long,  thin  cylinder.   The  perturbation  solution  for 
transverse  deflection  is  assumed  to  be  a  Fourier  cosine  series 


Z-*i       n 


w'  =  >    A  (t)cos  n6  (4.1) 

n=0 
4.2   Basic  Cylindrical  Shell  Equations 

4.2.1   Differential  Geometry 

Defining  x  as  n   and  6  as  n  ,  as  illustrated  in  Figure  4,  the 
1  2 

position  vector  equals 

—    /\        /N       /\ 

r  =   xi  +  asin6j  +  acos8k    .  (4.2a) 

The  principle  curvature  magnitudes  are 

* 

k   =  0  (4.2b) 

x 


and 


* 

Ka   =    1/a  (4.2c) 


and  the  differential  area  is 

ds2  =  dx2  +  a2d92  (4. 2d) 

where, 

A  -  1 
and  . 

B  =  a 
4.2.2   Strain  Curvature  Relations 

Equations  (2.7  a-f)  reduce  to 
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and 


and 


1  9 

e      =  u,      +  -  w,    z  (4.3a) 

x  x        2        x 

i  w'e2 

ee  =  I(v'e  +  w  +  -2r-)  (4-3b) 

u'o  w'xw'e 

y      =  +   v,       +  (4.3c) 

0  a  x  a 


k      =   -w,  (4.3d) 

X  XX 


Ke  =  \  (v'e  ~  w'ee)  (4'3e) 


2k    a  =  -   (v,      -    2w,    .)  (4.3f) 

x8        a  x  x6 


The   time   derivatives   of    (4.3a-f)    are 


e      =  u.      +  w,    w,  (4.4a) 

X  X  XX 


^e  =  I  (^e  +  "  +  Iw'e;'o)  (4'4b) 


•        i  •  •  l    *  • 

T      =  —  u,      +   v,      +  —   (w,   w,    +w,    w,    )  (4.4c) 

0        a        6  xa  x6        9x 


k      =   -w,  (4.4d) 

X  XX 


T  (^'e  "  ;'oe)  (4-4e) 

a 


2k        =  -   (v,      -   2w,n    )  (4.4f) 

x9        a  x  Gx 

4.2.3   Equilibrium  Equations 

Equations  (2.9  a-c)  yield  in-plane  equations 


\  ap+NQa+aN=0  (4.5a) 

x  x0,6  x,x 

and 
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apfl   +   aN  +  N+-M+M.        =0  (4.5b) 

6  x6,x  6,6        a      0,6  xQ,x 


and  transverse   equation 


ap      -   N.    +    a(N  w,         +  N        w,    ) 
z  0  xxx  x,xx 


+  I(V'ee  +  Ne,ew'e)  +  2NxOw'x( 


+  N  w,n   +   N  w,      +   aM 

xG,x      6  x8,6      x  x,xx 

M 

+  -^-  +   2M   fi       .   =   0  (4.5c) 

a  x0,x8 

4 . 3      Buckling   Solution 

For  the    axisymmetric    loading  of   an   infinitely   long   cylinder, 
the    following   assumptions   are   stipulated: 

(1)  loading  is    axisymmetric,    thus   the   dominant  solution  has 
no   8   dependence   and  p„   =   p     =0 

X  Ke 

(2)  plane  strain  conditions  exist,  N  =  -—- 

x    2 

(3)  transverse  displacement,  forces  and  moments  do  not  vary 
significantly  with  the  longitudinal,  x,  direction,  hence  both  the  dominant 
and  perturbation  solutions  are  independent  of  x. 

The  equilibrium  equations  simplify  to  dominant  equation  of 

p =  0  (4.6a) 

a    a 

and  perturbation  equation   of 


V'ee  +  M'e,ee    =  °  (4-6b) 

The  dominant  rates  of  membrane  strains  are 

F  =  0  (4.7a) 

x 
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efl  =  -  (4.7b) 

0        a 

and 

7=0  (4.7c) 

0 

and   the   rate   of   curvature   change  become 

•  •  • 

K   =  0       Kn=0       Kn=0 

x  0  x9 

Assuming  inextensible  buckling  the  perturbation  membrane  strains 

are  zero  and  the  curvature  changes  are  simplified  to 

k'   =  0  (4.7d) 

x 

K'e---^  (w«  -w',eQ)  ,    (4.7e) 

a 

K»    =  0  (4.7f) 

x8 

4.3.1      Linear  Vis co-Elastic   Case 

Substituting  the   constitutive   relations    (3.3    a-f)    into 

(4.6a)    reveals    a  dominant   transverse   displacement   rate   of 

2 
j_       a  p 

w  = .  (4.8) 

a 

The  perturbation  solution  (4.6b) ,  with  the  assumption  of 

inextensible  buckling  (e1   =  e*   =  E1  „  =  0)  ,  becomes 

x     0     x0 

aww'/ee  -  ?/a(w'/ee  +w',eQee)  =  0      (4.9) 

Substituting  the  assumed  solution  for  w1,  (4.1),  into  equation 

(4.9)  yields  a  first  order,  ordinary  differential  equation  of  constant 

coefficients . 

XA  +  (n2-l)A  =  0      (4.10a) 
n  n 

where 


-26- 


Thus: 


where 


A     =   amplification    function 
n 


•aaw 


A     =   C      expf — 2£ ] 

n  n  \(n2-l)    / 


(4.10b) 


(4.11) 


C  =  an  arbitrary  constant  dependent  on  the 
initial  imperfections. 


Defining 


♦  -  xt 

one   can  plot   A  /C      versus    t      to   discover  the   critical  mode    (n      . , ) . 
n      n  c  crit 

Figure   5   reveals    that   the   elliptical  mode,    n  equal   to  2 ,    is    the    critical 
mode.      The   growth   of   the    critical   amplification    function   is   displayed 
in  Figure   6 . 
4.3.2      Non-Linear  Vis co-Elastic  Case 

Combining  the   equilibrium  equation    (4.6   a-b)    and  the   non- 
linear  constitutive   equations    (3.8   a-f)    generates   a  dominant   transverse 

displacement   rate   of  ,      , , 

^_       a  p  a* 

(4.12) 


w  = 


a  p  cp 
g 


The  perturbation  equation  is  transformed  to 

fiww',ee  "  S/a<w',ee+w',e099)  =0      (4.13) 
Note,  if  creep  exponent  M  equals  one,  equation  (4.13)  reduces  to  the 
linear  equation  (4.8). 

Substituting  the  assumed  solution  for  w',  equation  (4.1),  into 
equation  (4.13)  yields  a  first  order  ordinary  differential  equation  of 
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constant  coefficients. 


where 


and 


hence 


XA      +     (n2+l)A      =    0  (4.14a) 

n  n 


A     =   amplification    function 
n 


X    =-^p-  (4.14b) 


(n2-l)/ 


A     =    C      exp  J I  (4.15) 

"  n  \(i- 

where  C  is  an  arbitrary  constant  dependent  on  initial  imperfections 
and 

T        =    Xt         . 

c 

As   in   the    linear   case,    the   elliptical   mode   is    critical. 

In   the    formulation   of  the    analysis,    it  was    assumed  that   the 
effective   stress    remained   constant,    however,    as    the  perturbation   solu- 
tion  of   transverse   displacement   grows,    the   effective   stress   increases. 
To  overcome   this  problem  one   can   divide    the   deformation  process   into  a 
finite  number  of  time   segments   in  which   the   effective   stress   is    considered 
constant.      This   operation   generates    a  piecewise    approximation   of   the 
transverse   deflection. 

Plane    strain   conditions    allows   the   effective    stress   to  be   ex- 
pressed  as 


o*   =   0.866    (a.   +   a'    )  (4.16) 

0  u 

where 
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p2a 

o„   =  -P-  (4.17a) 


which   reduces   to 


Defining 


a* 
0 


the  perturbation   solution  becomes 


4         (1-M)     • 


0.      _  2  Sla* <!-»)«.         .  ,4.17b) 

6  a2 


♦  - 


and 


c       (M-l) 
w'    =  C     exp    {-—  i>  )    cos26  (4.18a) 

n  3 


. (M-l)     T 
w'    =~ ~PW'  '  (4.18b) 


The   effective   stress  becomes 


(1±,.(M-1)  ,a.  /Tc.  (M-l)\    ] 

1+6^  (— )    a*expl^pjj  I    J 


p   a 
a*   =    .866  -—-[  H-6ip "  a*expl-^x"  I  (4.19) 


where 


a*   =   C  /a    (initial   out   of   roundness   pa- 
ri r 

rameter   <<    1) . 

The  initial  effective  stress  is 

x   =  X  t 
c     0 


X   -  initial  value  of  X 
0 

p  a 

0*   =  .866  -—■   (1  +  6  (^)  a*)  (4.20) 

0  h  h 
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The  duration  of  the  time  segment  is  dictated  by  the  growth  of 
effective  stress.   When  the  effective  stress  increases  by  a  given  per- 
cent, i.e.  ten  percent  (^  =  1.1) ,  a  new  time  segment  using  the  new  value 
of  a*  is  initiated. 

As  an  example,  taking  (a/h)  equal  to  10  and  a*  equal  to  .005, 
the  amplification  function's  accelerated  growth  causes  buckling  at  x 
equal  to  2.59,  as  graphically  portrayed  in  Figure  7.   Inspection  of 
Figure  7  shows  that  the  amplification  function  growth  is  very  slow.   It 
is  only  over  the  final  10  percent  of  the  total  lifetime  that  the  ampli- 
fication function  becomes  substantial,  indicating  catastrophic  buckling. 
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5.   COMPLETE  SPHERICAL  SHELL 

5. 1   Introduction 

The  succeeding  chapter  presents  the  creep  response  for 
the  complete  spherical  shell.   The  perturbation  solution  for  trans- 
verse deflection  is  assumed  to  be 


00  n 

y       (B   (t)P  (cos<}>)  +  J     E 


y     B   (t)pm(cos(j))co£ 
n=0  '  iSL 


where   P      and  P      are   Legendre   and  associated  Legendre  polynomials   of 

degree   n   and   order  m. 

5. 2      Basic   Spherical   Shell   Equations 

5.2.1      Differential   Geometry 

Stipulating   <j>   as   n      and   6   as   n    ,    as   exhibited  in   Figure   8, 

1  2 

defines    the   position  vector  equal   to 

r   =   R(sin<j>   cos9   i   +   cosA   sin8    j    -   sine})  k)  ,       (5.2a) 


the  principle   curvature   magnatudes   equal   to 

1 

<i>       '    R 


*  1 

K  a     =   TT  (5.2b) 


and 


6   '  R   ' 
and  the  differential  area  equal  to 


*  1 

K0    =  ■=■      ,  (5.2c) 


ds2   =   R2d(j)2   +   R2sin2Ad02  (5. 2d) 


where, 

A  =    R 
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and 

B   =    R  sine}) 
5.2.2      Strain   Curvature    Relations 

Equations    (2.7   a-f)    simplify   to  membrane   strains 

1    I  (W'<i>)2\ 

cn   -  —  [  v,„csc6  +   u   cotci   +  w 

e      r  \     8 


""e|2 

and 

Yo  =  r(U'9CSC*  +  v'*  "  V  COt^ 


csc2<H  (5.3b) 


and  curvature   changes 


and 


+  — | — -  csc2<j>)  (5.3c) 


k      =  (u,      -  w,       )  (5.3d) 

csc<\>     .                                         ,.            COttj)  .  ,               ,_     _    . 

Ke  ^  -T"  (v'o  "  w'eeCSC(j))  +  T~  ^V         (5'3e) 

R^  Rz 


„  csccj)    .  .  1     .  ,  , 

2k        = (u,       -  w,       )    +  (v,    -v   cot* 

<t>9  R2  e  $9  R2  <$> 


w,.,csc<J>   +   2w,  ncsc4>cotcf>)  .       (5.3f) 
t)(p  U 


The   time   derivatives    of    (5.3   a-f)    are 


^    =    £(%    +    "    +    R^V  (5'4a) 

*  1    *  •  •         esc    cb  • 

ee  =  p(v'0CSC(t)  +  u  cot$  +  w  +  — R — w,Gw,o)     (5-4b) 

•  1  /  •  •  •  esc   (!)    * 

Yq   =  -(esc*   u,Q   -   v/(j)   +  v   cot*   +  -r- (w^w,e 


+  wf<j)w,e))( 
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and 


K     =  —   (u,.    -  w       .)  (5.4d) 

<P         R2  9  99 

csccb    ,•              •                ,,  cot*      •    •      .  ,_     .    . 

Ke  =      o     (v'e  "  w'eecsccf,)  +  ~~T  (u"w'*)        (5*4e) 


•  csc<l)       *  *  *  * 

2k        =  (u,      -  w,      )    +  —   (v,      -   v  cot? 

9°  R2  9  99  R2  9 


-  w, ..cscQ   +   2w, ^cscOcoto)  (5.4f) 

99  0 


5.2.3      Equilibrium  Equations 

Equations    (2.9    a-c)    simplify   into  in-plane   equations 


Rpd)   -   N0cot9+  11    n    ncsc(j)  +   Nx     ,    +   N,cot9 

96f  8  9,99 

+  —    (M,     .    +   M,cot9   -   M.cotcj)   +   M,_    ncsc9)    =   0  (5.5a) 

R9,9        9  0  90, er 


^9   +   N90,9   +   2COt*N90   +   CScn,0 

+  ±    (csc9Mo>0   +   M^^   +   2M<{)0cot9)    =   0  (5.5b) 


and  transverse   equation 


R2p2   +  ^(w^cot*   +   w,w)    +   NeW'09CS^( 


R(N9  +  No)  +  2csc*N90w'69  +  \,9W'9 


+  csc29Ne^0w,e  +  csc9(N     wrQ  +  N^w,^) 


+   M9,99  +   2COt*M9,9   "   %   +   CSC2<|,M0,09 


+    2M.„    ncot9csc9   +   2Mjn     .  „csc<{j    -   ML     .cot9 
90,0  90,90  0 1 9 


+   M      =    0  (5.5c) 
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5. 3      Buckling  Solution 

In   the    analysis   of   the   complete   sphere,    it   is    assumed   the 

dominant   solution  has  point   symmetry,    loaded  in   a  uniform  pressure 

field.      As   a   consequence,    there   is   no  dominant    0   or   (J   dependence, 

P6   =  P<j>  =   0,    M,  .  _   =   0,    M,    =   M„ ,    and  N„   =  N,.      The    dominant   transverse 
90  9  8  0  9 


equilibrium  equation   simplifies   to 


Rp      -    2N.    =   0  (5.6a) 

rz  <j> 


The   in-plane  perturbation   equilibrium  equations    reduce   to 


M'        .    +    cotcfuM'       -    cot<J>M*       +    csc<j>M'         a  =   0  (5.6b) 

9,  <p  q>  6  (j>0 ,  0 


and 


csc<J>M'       n    +    M±n     .    +    2cot<j>M*         =    0  (5.6c) 

0,0  90,9  90 

The   transverse  perturbation  equilibrium  equation  becomes 

VW''4>*   +    COt*W''9   +    CSc2<J,W''O0)    +   M9,99 

+    2cot9M'  -    M'       +    csc2<f>M'       nr>    +    2cot9csccf>M '     .   . 

9,9  9  T       8,60  90,0 

+    2csctj)M'     .     ..    -    cot<J>M'        .    +  M»      =  0  (5.6d) 

90,90  0,9  0 

The  dominant  rates  of  membrane  strains  and  curvatures 

changes  are 

•     »     —  • 

79  =  F0  =  I        '  F90  =  °     (5-7a"C) 

and 

•  •  • 

^  -  7Q   =  0        ,  iT+e  -  0     (5.7d-f) 

Assuming  inextensible  buckling  the  perturbation  rates  of  membrane  strains 
are  zero  and  the  rates  of  curvature  changes  can  be  expressed  as 
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and 


K'      =  -  —   (w'    +  w\  ..)  (5.8a) 

<P  R2  M 

k'      = (w'    +    csc2cj)w'       Q   +    cot^w',.)  (5.8b) 

V  2  oo  9 

2 

2k'     .    =  —    (cscicotiw1 , n   -    cscow' , . A)  (5.8c) 

<J)9         r2  6  (j>e 

If   one   examines   the  perturbation   expressions   derived  for   cur- 
vature  changes  when   inextensible  buckling   is    assumed    (equations   5.8a  - 
5.8c),    one   discovers   that   they    are   identical   to   the    curvature    change 
expressions   derived  by  Jones    and  Ann    (1974).      By   comparing  the    formu- 
lation  of   the   transverse   equilibrium  equations   of  Walters    and  Jones    and 

W 
Ahn,    it  is    found  that  Walters    and  Jones    assumed  the  —  term  to  be   neg- 

R2 

ligible   in   the    curvature    change   expressions,   while   Jones    and  Ahn   include 

it.      The    curvature   change   expressions    are    assumed   to  be    a   function   of 

u,    v,    and  w  by  Walters    and  Jones,    while   Jones    and  Ahn   assume   curvature 

change   expressions    are    just   a   function   of  w.      If,    however,    inextensible 

buckling  is    assumed,    one    can  express   the   in-plane   displacements,    u   and 

v  as    functions   of   the    transverse   displacement  w   only,    and   arrive    at 

the   expressions    derived  by  Jones    and  Ahn.      Thus,    it  is    reasonable   to 

assume   that   the   transverse   equilibrium  equation   of  Walters    and  Jones 

must  be   modified  by   inclusion   of  the   term 


(5.9) 


which   is    generated  by   the  —  term. 

R2 
Hence,    the   modified  transverse    perturbation  equation  becomes 


NA(W'»A,     +    cotcjw',,     +    csc2<j)w'  ,       )     +    M 
<j>  <p<j>  <p  80  (, 
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+    2cot<pM'        ,    +    2M'       +    csc2cf)M'       nn 

9,9  0  e,00 

+    2cot9csc9M'  ,  .    .    +   2csc<f>M'  ,_.    -    cotcpM'       ,    =  0  (5.10a) 

(p0  ,  0  (p0  ,  (p0  0,9 

The   in-plane  perturbation   equations,     (5.6  b-c)    are   rewritten 
in   terms   of  w   as 

f"''«p   +  W''W   "  |cso29cot(}>w-/eQ   +  fcsc2^',^   =   0  (5.10b) 

and 

5  o      — 

cotcpw' ,        +    (—  -   4cot^4>)w, 

o    •  5 

+   csc^fpw',  +  —  csccpcotcpw,  =   0  (5.10c) 

To   remain   consistent,    the  preceding   alteration   is   only   valid 
if  equations    (5.10  b-c)    are   either  equal   to   zero   or  nearly   zero.      The 
equations    can  be   nearly   zero   since   the   transverse   equation    (5.10a)    is 
considered  the   more    dominant  equation. 
5.3.1      Linear  Visco-elastic   Case 

The   dominant   transverse   displacement   rate   is    achieved  by   sub- 
stituting  the    linear   constitutive   equation    (3.3a)    into  the   equilibrium 

equation    (5.6a) 

2 
_1         a   P 
W    =    -r-^  (5.11) 

3a 

Utilizing   the    following   Legendre    identities 

Pn,  ,  .    +    cotcpPn,.    +    (csc2<p   +    An)Pn   -   0  (5.12a) 


and 


where 


Pmn,JA   +    cot(pPmn,A    +    (m2csc29    +    An)  Pnn   =   0  (5.12b) 

99  9 


An  =  n(n  +   1) 
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The  perturbation   equation    (5.10a)    can  be   recast  to 


?2,7i     _     /n2     ,    i  /n2 


$Vzw'    -    (V^    +  -)  (V-    +    2)w'    -    0  (5.13a) 


where 


V2(    )    =    (    ),        +   cot(j)(    ),      +   csc2(    ),ee  (5.13b) 

and 

.  3aaw 

$  =   _  

2? 

Substituting  the    assumed  solution    for  the  perturbation   trans- 
verse  displacement    (5.1)    uncouples   the    amplification    function   in   the 
form  of   ordinary   differential  equations    of   constant   coefficients 

$AnBrt      +    (t  -    An) (2    -    An)B_      =    0  (5.14a) 

On  2  On 

$AnB        +    {—  -    An)  (2    -    An)B        =   0  (5.14b) 

mn  2  mn 

Equation    (5.14b)    is   identical  to  equation    (5.14a)    if  m  is   understood  to 
include  the   case   m  equals    zero. 

Thus,    providing  p      is    constant 

I  OAnt  (  ._    ,_, 

B        =   D        exp     ~ )  (5.15) 

mn  inn  J    ,3 

(     2    ~  "  1 

where   D        is    an   arbitrary   constant   dependent   on   initial   imperfections. 

When    x      is   defined   as    4>t   a  non-dimensional  plot,    Figure   9,    can 
be   constructed   to  exhibit   that   the    critical  mode   is   n  equal  to  two. 

If   the    function 

2 

w'        X~^     m 

- —  =    >       P       (cos<j))cosme  (5.16) 

is  calculated  for  a  number  of  points  of  various  <■>   and  6  values  using 
computer  program  A  of  Appendix  (A) ,  the  two  dimensional  contour  surface 
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of   Figure    10   is    generated.      On   the    contour  map,   positive  values   in- 
dicate   an   outward  elevation   of  the   mid-surface    from  the   mean   radius,    and 
a  negative   value    indicates    an   inward  depression   of  the   mid-surface    from 
the  mean   radius,    due   to   the  perturbation   solution.      Since   the   elliptical 
mode  has   been  proven   to  be   the   most   critical  mode,    it  will  be  henceforth 
assumed  that  the   initial   imperfections    are   of  the   elliptical   shape. 
5.3.2      Non-linear  Visco-elastic   Case 

Substituting  the    constitutive   equation    (3.8a)    into    (5.6a) 
produces    a  dominant   transverse   displacement   rate   of 

w  =  — ^ '  <5-«> 

In  a   method  similar  to  the  linear  case  (eqn.  5.12a) ,  the 
perturbation  transverse  equilibrium  equation  is  transformed  to 

AV2w'  -  (V2  +  J)  (V2  +  2)w'  -  0     (5.18a) 
with 

A--f^  ■    (5.18b) 

When  the  curvature  change  equations  (5.4  d-f)  are  combined  with  consti- 
tutive equations  (3.8  d-f)  and  then  substituted  into  equation  (5.10). 

Using  the  proper  Legendre  identities  (5.13  a-b) ,  and  assumed 
solution  for  w1  (5.1),  the  uncoupled  expressions  for  the  amplification 
function  become 

AAnB    +  (|-  -  An)  (2  -  An)B    =  0     (5.19a) 
on     2  on 


and 


AAnB    +  (|-  -  An)  (2  -  An)B    -  0     (5.19b) 
mn    2  mn 
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Thus,  providing  p   is  constant 


where 


Anx 

B    =  D   exp|  — |  (5.20) 

(f  -  An)  (2  -  An) 


D   is  an  arbitrary  constant  dependent  on 

mn 


initial  imperfections, 


and 


x      =   At 
s 

The   critical  mode,    as   in   the    linear  example,    is   n   equal   to  two. 

As   with   the   infinitely    long   cylinder,    the    assumption  that   the 
effective   stress   remains    constant   must  be   modified.      Effective   stress   is 
assumed  to  be  piecewise   constant.      Time   is   divided  into   a   finite   number 
of  segments   within  which   the   effective   stress    is    assumed  to   remain   con- 
stant.     When   the   effective   stress   increases  by   a  given  percentage,    the 
segment  is   terminated,    and  a  new   segment  with    an   increased  value   of 
effective   stress    is    initiated. 

From  equation    (3.5b)    the   effective   stress   equals 

a*  =  -L  ((a      -  aQ)2   +  a.2   +  a   2  +  6x.fi    J  (5.21a) 

^2    V      ♦  6  *  6  (P9/ 

where,    from  equations    (1.2    a-f   and   3.8   a-f) ,    the   stresses    are   expressed   as 


o.    =   o.    +   o    '  (5.21b) 

$  $  <j> 


ae  =  °e  +  V  (5-21c) 


Td>e  =  T'*e  (5'21d) 
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From  equations  (5.7  a-f)  and  5.8  a-c)  expressions  (5.21b  -  5.21d)  can 

be  expressed  as 

_    _    P  R 

o,  =  o   =  -~-  (5.22a) 

(p  6    2H 

where  H  is  the  shell  thickness  of  the  sphere, 

4  (1-M)     •  1    • 

cr'      =  ^-Kho*^        (e'      +  ±  e»    )  (5.22b) 

q>  3  <p  2  o 

rewritten   in   terms    of  w 

2  KH      +(1-M)  .3       ,  •  W'  '(f) 

+  — j2^-  csc2<f>)     (5.22c) 
which   can  be    further   simplified   using  Legendre   identities    (5.12   a-b)    to 

KH     (1-M) 
a'   =  — —  o*K        '   (3w-  -  w'   .)  (5.22d) 

*    3R2  ** 

and  in  a  similar  manner 

K  II       .  ( 1-M)  ,  3 •  • .  •  .  .  •  .  2  a 

o'Q   = a*  (jw1    +  w,    cotiji   +  w'  ,Q6csc^ 

R2  9 

+  T"''aa'     (5.22e) 
2  <J>q> 


which   is   simplified  to 


KH    (1-M) 
0'   =  -^-  a*VX  iU  (9w'  +  w,  .  ,)  (5.22f) 

e  •  3r2  M 


Finally  the  shear  stress  is  written  as 


K    (1-M)  • 


which,    by   substituting  equation    (5.8c),    is 


2KHa*(1-M)      • 

T1         =   (w1  ,  ,cscJ)cot(J)    -   w',,_csc(f))  (5.2     l) 

4>0  2  ♦  <P0 
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Since  the  perturbated  strains  depend  only  on  curvature  change 
rates  (5.8  a-e) ,  which  are  all  expressed  as  functions  of  w1,  it  is 
postulated  that  the  spatial  point  on  the  sphere  which  exhibits  the  maxi- 
mum perturbated  transverse  deflection  rate  will  be  the  point  of  maximum 
strains,  thus  stresses-  To  locate  the  point  of  maximum  curvature  rate  (w') 
one  can  utilize  a  precedure  to  study  maxima  and  minima  of  functions  of 
several  variables  outlined  by  Hilderbrand  (1962). 

(a)  Evaluate  w1,   and  w1,   and  then  set  equal  to  zero. 

<p  0 

(b)  Evaluate  w',^      ,      w\Qe    and  w',^ 

(i)       if  w'     .w'    „   >  w',./      ,    and  w'     ,    >    0      ,    a   relative 
99      06  98  99 

minimum  exists, 

•        •  •  • 

(ii)   ifw,..w'_>w'_   ,  and  w1 , ,  <  0   ,  a  relative 

99   60      96  99 

maximum  exists. 
From  equation  (5.1)  the  critical  value  of  w' (n  =  2)  is 

w*  =  B   (t) (p  (cos9)  +  COS0P1  (cos9) 

9m     \  q  o 


2m 


+  cos20P2  (cos(j)))  (5.23a) 

9      / 


with 


Thus , 


P  (cos9)  =  -  (3cos29  +  1)  (5.23b) 

2  2 

dP2 (cos9) 

P1  (cosa)  =  sind)  — — — —  =  3sin9cos9  (5.23c) 

2  d(costj)) 

d2P2  (coscj)) 

P2  (cos9)  =  sin29  -  3sin29      .  (5.23d) 

2  d(cos9)2 


•  9  9 

w',9   =   -3cos<|)sin9   +   3cos6(cos^9   -   sin^9) 

+   6cos20sin9cos9  (5.23e) 
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*  7 

w',   =  -3sin8sin(j)cos<})  -  6sin26sin  9  (5.23f) 

8 

Both  of  these  expressions  (5.23  e-f)  are  zero  when  9  and  G  equal 

7T/2. 

Evaluating   the   second  derivatives 

•  0  7 

w' ,  .  .    =    3  (sin   9   -   cos    A)    -    12cos8sin9cos9 

0  9 

+   6cos28(cos    9   -    sin    9)  (5.23g) 

7 

w'    =  -3cos8sin(j)cos9  -  12cos2Qsin  9  (5.23h) 

Uo 

and 

*  9  9  1 

w1,        =   +3sin8(cos    9   -   sin    9)    -    12sin2  9cos(J>3inrj>    (5.23i) 
98 

at   9   =   tt/2    and   9   =   tt/2   one    finds 

w\         -   9  (5.24a) 

<P<P 

w',0Q   =    12  (5.24b) 

w',.fl   =    -3  (5.24c) 

90 

•  •  • 

Thus,    w',,.w',A.    >   w',.„    and   a   relative   minimum  exists    at   9 
99         88  98 

equal   to  tt/2    and   8   equal   to  tt/2.      Noting,  however,    that   the   deflections 

rates   herein    are   of   a  negative   sense    (B        is   always   negative) ,    the   point 

2m 

9  =   tt/2,    8   =   tt/2    is    actually   a   relative   maximum.      Further,    this    conclusion 
is  in  agreement  with    the   maximum  of  the   deflection    field   illustrated  in 
Figure    10,    which  was    calculated  numerically. 

Substituting  values    (5.24   a-c)    into  equations    (5.22   d,    f,    h) 

yields 

a.      .         1*L0*U-M>£       (t)  (5.25a) 

?  l    R2  2m 
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7  KH      ^(1-M);       ... 
cr»      =  —  —  a*  B      (t) 

0         2    R2  2m 

T..„    =    2KHa*U-M>£ 


<t>G 


Rz 


2m 


(5.25b) 


(5.25c) 


From  equation    (5.15) 


2    T  T 

B'       (t)    =  — -f-exp  -f- 
2m  3      t  3 


3   ^    .R,  3   P2       .  (M-l)  ,Ts. 


the  perturbated  stresses   simplify   to 


(5.25d) 


a' 


3  .R.  o   .    2      .  s. 

—  (— )ZRX  -—  exp(— -) 

4  H      k    *   3 


(5.25e) 


P  T 

a'  =  21  ,R.  ?„*  «     /  s,    _  . 


(5.25f) 


T'*6=  3(f)R*^exP(i)  =  -40- 


where 

R*  = 
Hence  equation  (5.21a)  becomes 


(5.25g) 


maximum  deviation    from  radius    of  imperfections 

radius 


a*  =  —  |(8a\)2  +    (a 
/2    L         ♦ 


+  a  '  J2   +    (o,   -   la')2 


zh 


+   6(-4a    )z         (5.26a) 


simplified   to 


a*  -  —     2a   2   -   120,0'  ,   +  200a1    2] 
./T  I  $      *  *   J 


Substituting  expressions    (5.22a   and   5.25e),    the   effective 


(5.26b) 


stress   equals 


1    PR 

a  H 


1        9    ,R.  „.  ,    s, 

2"2    (i)R*  ^^ 


+    112.5  (£)2R*2exp2(-f) 
H  3    J 

The  initial  effective  stress  (t  =0)  equals 

s 


(5.26c) 
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Defining 


0  v/J 


4>  = 


H 


M;^""0V 


(5.26d) 


(5.27a) 


and 


x      -At 
s  0 


A      =    initial   value    of   A    (eqn.    5.18b) 
0 


(5.27b) 


the   effective   stress    can  be   recast   as 


a*  = 


/2    H 


1  9     ,RX      ..M-l  sy 

2  "    2     (H}  R **         eXP  ~3 


M-l 


+    112(f)2R.2^«-D2cxp2^! 


M-l 


(5.27c) 


As    an   example,    with   R/H,    R*,    M,    ijj     equal   to   10,    .005,    3, 

0 

1.1   respectively,    the    growth   of   the   amplification    function   is    as    shown 
in   Figure    11.      Note   that   the    amplification    function   stays    relatively   small 
until   the    last  ten  percent   of  the    lifetime   of   the   sphere.      As   expected, 
the   time   segment   duration   decreases    as   time   increases. 

It   should  be    further  noted   that   the   in-plane   equilibrium  equations 
(5.10  b   &   c)    are   equal  to  zero   at   the  point   of  maximum  stress    and  trans- 
verse   deflection    (<J)   =    it/2  ,    0   =    ir/2)  . 
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6.   DISCUSSION 

6. 1   Infinitely  Long  Cylinder 
6.1.1   Comparison  to  Past  Results 

The  most  convenient  basis  of  comparison  of  the  analysis  herein 
with  past  investigations  is  the  predicted  critical  time  of  creep  buckling. 
The  time  behavior  of  the  initial  imperfection  growth  is  such  that  it  will 
become  large  and  grow  rapidly  only  for  time  near  to  the  critical  time.   As 
an  aid  for  the  reader,  Appendix  B  has  been  provided  to  transform  the  nota- 
tion of  past  analyses  into  the  notation  used  herein. 

Hoff  et  al  (1959)  predict  a  critical  time  of 

t    =  ^ (6.1) 

k   h 

It  is  at  this  critical  time  at  which  the  amplification  of  the 
initial  imperfections  becomes  indefinitely  large.   The  herein  results  have 
been  compared  to  Hoff  in  three  ways  by  varying  the  values  of  a/h ,   a*,  and 
M  in  Figures  12,  13,  14  respectively.   The  comparison  reveals  favorable 
agreement  for  smaller  a/h  ratios  (a/h  <  25)  with  a*  and  M  equal  to  .005 
and  3  respectively;  and  small  imperfections  (a*  <  .003)  with  a/h  equal  to 
10  and  M  equal  to  3.   For  large  a/h  ratios,  the  Hoff  estimate  of  buckling 
time  is  approximately  an  order  of  magnitude  greater  than  the  presented 
analysis.   The  Hoff  solution  is  less  sensitive  to  larger,  gross  initial 
imperfections,  predicting  a  collapse  time  of  five  times  greater  than  the 
herein  analysis  when  gross  imperfections  are  assumed. 

The  plot  of  varying  creep  exponent,  M,  (Figure  14)  shows  excellent 
agreement  for  M  equal  to  one  and  seven  and  an  estimate  five  times  less  than 
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Hoff   for  M  equal   to   three    and   five.      In   Figure    14,    a*   and   a/h  equal 
.005    and   10    respectively. 

The    actual   time  history   of  the   two   amplification   curves    are 
nearly   identical,    as   exhibited  in   Figure    15.      The    amplification   of   im- 
perfections   remains    very   small   until  the    last   ten  percent   of   the    life- 
time  of   the    cylinder. 

The    discrepancy   in  predicted   critical   time   of  the   two   analyses 
do  not   seem  too  serious   in   the    light   of   the    fact   that   creep  buckling 
times    are   of  the   order  of   105   hours,    so  that  excellent   agreement  is   ob- 
served  for  the    first   10^   hours. 

Berman  et   al    (1974)    presented  the  particular   case   of  stainless 
steel   304,    which  has   the   material  properties    at   1050   F   listed  in   Ap- 
pendix C.      One    finds    a  great   overestimation   of  Berman  when   compared  to 
both   the   analyses   of  Hoff   and  herein.      Berman  predicts    a  buckling   time 
for   cylinder   of   a/h   equal  to   10    and   a*   equal   to    .05   of   1.743   x   102t*p   "5-G789 
while  Hoff  predicts    3.6852    x   102°p      ^"6789^    and  the   herein   investigation 
estimates    1.826   x   102°p  "5-6  789^      Both   the   Hoff   and  herein   analysis 

2 

predict   a  buckling  time   of   approximately   four  orders   of  magnitude   sooner 
than   Berman.      The   author   feels   that   the   calculated  time    of  Berman   is   in 
possible   error,    since   he   uses    an   out-of-roundness    ratio   of    .05,    a  very 
gross   out-of-roundness    for   any   realistic  engineering  design.      Further 
Berman    fails    to  update   the   effective   stress  with   time.      The  herein 
analysis   does   predict   a  buckling  time   of   8.3   x   1022p      5-6789    ^f   a*   ^s 
.005,     (an   order  of  magnitude    less)     and  the   effective   stress   is   not   up- 
dated with   increasing   time. 
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The   numerical    analysis    of  Samuelson    (1970)    deals  with   the   creep 
response   of   finite    length   cylinders   subject   to  both   non-uniform  and 
uniform  loads.      Thus,    higher  harmonics    come   into   an   influencial   role   in 
the   creep  process.      Further  the   transverse   deflection   is   dependent   on 
the    longitudinal,    x,    coordinate 

CO 

V    =    /  ,    A    (t)cosnGsin  —  (6.2) 

*—i>      n  L 

n=0 

Samuelson  predicts    a  buckling  time   of   approximately   100 
seconds    for   a   cylinder  with    a/h,    L/a  equal  to  27   and  2    respectively. 

A  similar  buckling   time   is    found   for  the   same    cylinder  by   the  herein 

K 
analysis   if  M  equals    3,    a  nominal    a*     .002,    is    assumed  and  — *-     is 

p  ° 
a 

o 

assumed  to  be  of  the  order  10°.   Thus,  it  is  possible  to  state  that  the 
herein  analysis  predicts  a  buckling  time  of  the  same  order  of  magnitude 
as  Samuelson. 
6.1.2   Elastic  Response 

The  elastic  response  of  the  cylinder  can  be  obtained  by  sub- 
stituting the  elastic  part  of  the  constitutive  relation  (3.2a)  into  the 
transverse  equilibrium  equation  (4.6a)  to  yield 

a2pz(l-u2) 

w  =  — (6.3a) 

e      Eh 

The   elastic  transverse   deflection    f rom  Roark  (1965)    is 

a2P 

w  =  — -    (1   -   £)  (6.3b) 

e  _,        ,        Eh  2 

Roark 

If  Poisson's    ratio,    u,    equals    .3,    the    ratio  of  the  herein 

elastic  displacement   to  that  of  Roark  is    1.07. 
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6.1.3      Parametric  Study 

Due   to  the   dependence   of   the   tiine  history   of  the   transverse 
deflection   on   a  multitude   of  independent   variables,    radius   to  thickness 
ratio    (a/h) ,    out-of-roundness    ratio    (a*) ,    material   constants    (M,   k)    and 
loading    (p   ) ;    an   indefinitely    large   number  of  parametric   curves    are 
required  to   adequately   explore   the    creep   response   of   cylinders   under 
external  pressure.      This    task   is    greatly   simplified  by   the    computing 
of   transverse   displacement    (4.18a)    using   computer  program  B   of  Ap- 
pendix A. 

Figures    16,    17,    and   18   are  plots   of  non-dimensional   critical 

(CI:     2         | 
I,    versus    a/h    for   creep   exponent   M  equal  to   3,    5,    7    re- 
spectively.     In   each   case    as    a/h   ratio   for   a  given  a*   decreases,    the 
non-dimensional   critical   time    drastically   increases.      Further,    as    a* 
increases    for   a  given   a/h   ratio,    the   non-dimensional   critical   time   de- 
creases.     Of   the    two  parameters,    a/h   and   a*,    a/h   is   the   more   dominant 
variable.      That   is,    critical   time   is   more   sensitive   to   a/h   than   a*.      The 
same   conclusion  was    reached  by  Hoff    (1959) . 

Figure    19    illustrates   the   effect   of  varying  the   initiation 

point   of   a  new  time   segment   of   constant   effective   stress.      If   for  example, 

t  p  3 

CX  55 

M,  a/h,  a*  equal  3,  10,  .005  respectively,  the  value  ^ equals 

3.5  x  10  °    for  time  segments  which  are  initiated  when  the  effective 
stress  (4.19)  increases  by  fifty  percent  (ijj  =  1.5)  compared  to  a  value 
of  2.2  x  10  &  when  the  time  segments  are  initiated  when  the  effective 
stress  increases  by  five  percent  (ip  =  1.05),  a  ratio  of  1.6  difference. 

From  a  practical  ocean  engineering  point  of  view,  a  comparison 
of  the  two  most  common  creep  prone  materials,  aluminum  and  titanium,  used 
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in   the   ocean  is   in   order.      The   alloys   used   in   the    comparison   are   Aluminum 
6061-T6    and  Titanium  6Al-2cb-lTa-.08Mo.      The    room  temperature    creep 
constants   of  these   materials    are    listed  in   Appendix  C.      The   titanium  con- 
stants  where  experimentally    found  by   Chu    (19  72) ,   while   the    aluminum  con- 
stants by   the    author. 

The    first   comparison   is   non-dimensional   critical  time   verses 
the   initial  elastic  hoop   stress,    expressed   as    a  percentage   of  yield  stress 
The  yield  stresses   of  Titanium  621/. 08   and  Aluminum  6061-T6   are    120    and 
40  ksi   respectively.      The   two  materials    are    compared  with   two   a/h   ratios, 
10    and  50;    and  two  a*   ratios,    .01,    and    .001.      A  number   of   conclusions 
can  be   drawn.      The    lifetime    increases    for  both   materials,    approximately 
five   times    for   a/h   equal   to   10    and   forty-five   times    for   a/h   equal   to  50, 
with   an   order   of  magnitude   decrease   in   initial   out-of-roundness    ratio 
(a*).      Hence,    it   appears   to  be   more   important   that    cylinders   of   large 
a/h   ratios   to  be    constructed  to   closer  tolerances   than   cylinders   of  a 
smaller   a/h   ratio.      Titanium  is   the  better  performer  at   low   stress    levels, 
but   is   surpassed  by   aluminum  at  higher  stress    levels.      Imperfections  play 
a  more   important   role    at  higher  stress    levels.      Finally,    for   a  stress 
level    less   than    forty  percent   of  yield,    there   is    a  dramatic   increase   in 
lifetime. 

The   influence   the    varying   a/h   at   a  given   stress    level,   sixty 
percent   of  yield,    has    on   the   non-dimensional   critical   time   is   portrayed 
in   Figure   21.      Again,    a   crossover  occurs.      For  small   a/h   ratios    (a/h    <   20) 
titanium  is   the   optimum  material,    while    for   larger   a/h   ratios    aluminum 
exceeds   titanium. 
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Taking  another  tack,  if  a*  is  varyed  for  a  constant  a/h  (10) 
and  a  constant  stress  level,  fifty  percent  of  yield,  Figure  22  is  gen- 
erated, a  crossover  of  performance  is  observed  with  titanium  being 
surpassed  by  aluminum  for  grosser  tolerances  (a*  >  .007) . 

The  most  pertinent  collation  is  the  comparison  of  the  two 
materials  at  equal  pressure  levels.   As  an  example,  with  p  ,  a*,  a/h 
equal  to  2000  psi,  .001  and  10  respectively  titanium  out  performs 
aluminum  by  four  orders  of  magnitude,  (Figure  23).   Titanium,  with  its 
superior  strength,  is  able  to  operate  at  a  much  lower  percentage  of 
yield  stress. 
6. 2   Complete  Sphere 

6.2.1  Elastic  Response 

By   substituting   the   elastic  part  of   the    constitutive   equation 
(3.2a)    into   the   transverse   equilibrium  equation    (5.6a),    an   expression 
for  the  elastic   response    is    developed 

a2p2(l  -   u2) 

w      =  — : —  (6.4a) 

e  jEh 

The  elastic  deflection  given  by  Roark  (1965)  is 

a2p  (1  -  u) 

-  — fih (6-4b) 

Roark 

A  ratio  of  the  elastic  deflection  to  that  of  Roark  is  .86  of  u  is 

assumed  to  be  .3. 

6.2.2  Parametric  Study 

As  with  the  infinitely  long  cylinder,  the  dependence  if  the 
time  history  of  the  transverse  deflection  on  a  large  number  of  indepen- 
dent variables  requires  the  generation  of  a  large  assemblage  of  para- 
metric curves  to  adequately  investigate  creep  response  of  the  sphere. 
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Again,    the   task   is   simplified  by  numerically   computing  equation    (5.1) 
using   computer  program  C   of  Appendix  A. 

Figures    24,    25,    and  26    are    graphs    of  non-dimensional   time 

(  PA 

It  I  versus    R/H    for   various    creep  exponents,     (M  equal  to   3,    5,    7)  . 

V  c  K  / 

The   non-dimensional   critical   time   decreases    rapidly  with  both   increas- 
ing  R/H   and   R*,   being  most  sensitive   to   R/H. 

The  non-dimensional   critical   time   of   the   sphere   is    fairly 
non-sensitive   to  the    duration   of  the   time   segments,    as    illustrated   in 
Figure   27.       If  the   updating   of  effective   stress  (5.26a)   is    carried  out 

when   effective   stress   is    increased  by    fifty  percent    (ijj  =   1.5)    the 

t     P    3 
value   of is    1.6   x   10      .      If,    however,    the   effective   is    updated 

IN. 


t       P    3 

era 


when  it  has  increased  by  five  percent,  (^  =  1.05) equals  1.05 

K 

—  c 

x  10      ,    an   increase    of    fifty  percent. 

When   the   non-dimensional   critical   time    is    calculated   for 
various    initial   elastic   loadings,    expressed   as  percentage   of  yield 
stress,    for   the   two   alloys,    Aluminum  6061-T6    and  Titanium  621/. 08, 
Figure   28   is   manifested.      Initial   imperfections    are   quite   influencial 
on   the    critical   time,    causing   a   two  order   of  magnitude   increase    for 
a/h   equal   to   10    and  six  orders   of  magnitude    increase    for   R/H   equal 
to  50  with    an   order  of  magnitude   decrease    in   initial   imperfections. 
Thus,    revealing  that    for  spheres    of   large    R/H   ratios,    close   tolerances 
are   a  prerequisite.  Except    for  the    case   of  high   R/H    ratio   and   gross    im- 
perfections,   titanium  is   the   superior  material,    producing   a   longer 
lifetime.       Rapid  increase    in    lifetime    is    achieved   if  the    stress    level 
of   the   sphere   is    less   than   sixty  percent   of  the   yield  stress. 
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When  R/H  is  varied  for  a  given  stress  level,  sixty  percent  of 
yield,  titanium,  regardless  of  the  R/H  ratio,  out-performs  aluminum. 
This  superiority  is  even  more  obvious ,  as  demonstrated  in  Figure  29 , 
for  substantial  R/H  ratios  (R/H  >  30)  . 

If  for  a  constant  R/H  ratio  and  a  loading  of  fifty  percent 
of  yield,  the  magnitude  of  initial  imperfections  is  allowed  to  veer, 
Figure  30  is  developed.   Again,  titanium  is  always  the  superior  per- 
former, especially  with  small  initial  imperfections.   That  is,  one 
has  much  more  to  gain  by  machining  a  titanium  sphere  to  exacting 
tolerances  than  an  aluminum  sphere. 

At  the  same  pressure  level,  the  titanium  is  vastly  superior 
than  the  aluminum,  having  a  lifetime  of  approximately  7  orders  of 
magnitude  longer  with  p   equaling  2000  psi  and  R*  equaling  .001.   The 
obvious  reason  being  that  a  titanium  sphere  will  operate  at  one-fourth 
the  percentage  of  its  yield  stress  at  a  given  ambient  pressure. 
6.3   Comparison  of  the  Creep  Performance  of  the  Sphere  and  Cylinder 

If  the  critical  time  of  a  cylinder  and  a  sphere  of  the  same 
material  at  equal  pressure  levels  are  compared  by  plotting  the  ratio  of 
critical  time  of  the  sphere  to  the  critical  time  of  the  cylinder  for 
various  a/h  ratios,  initial  imperfections,  and  creep  exponents,  (Figures  31, 
32,  and  33)  it  is  discovered  that  regardless  of  the  initial  conditions, 
the  cylinder  has  the  inferior  creep  resistant  properties.   This  point  is 
especially  manifested  in  the  case  of  titanium,  M  equals  7.   The  oscil- 
lations observed  in  Figures  31,  32,  and  33  result  from  the  variation  of 
the  ratio  of  the  effective  stresses  with  various  a/la  and  a*  combinations. 
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The  sphere  is  more  sensitive  to  variations  in  initial  im- 
perfections, which  is  to  be  expected  since  the  cylinder  is  a  develop- 
able surface  and  is  less  constrained  than  the  sphere. 

The  sphere  has  excellent  creep  characteristics  for  stress 
levels  less  than  sixty  percent  of  yield  stress.   The  cylinder,  on  the 
other  hand,  the  stress  level  must  be  less  than  forty  percent  of  yield 
stress  for  favorable  creep  properties. 

Finally  titanium,  because  of  its  high  yield  stress  is  in 
every  instance  the  best  material  to  be  used  in  ocean  engineering  ap- 
plications of  both  the  sphere  and  the  cylinder. 
6. 4   Design  Criterion 

The  extensive  analysis  present  suggests  an  assemblage  of 
simple  design  guidelines  to  assist  the  engineer  in  developing  a  cost- 
effective,  creep  resistant  design.   In  most  ocean  engineering  shell 
designs,  the  major  constraint  in  design  is  considered  to  be  the  minimi- 
zation of  the  weight  to  displacement  ratio  of  the  structure,  while 
maintaining  a  maximum  stress  level  of  approximately  fifty  to  sixty  per- 
cent of  yield.   It  is  this  constraint  which  dictates  the  use  of  light- 
weight, but  creep-prone  materials  in  the  first  place.   Fortunately, 
this  criterion  is  satisfied  if  titanium  is  used  rather  than  aluminum  in 
pressure  vessels,  even  though  titanium  is  1.5  times  as  dense  as  aluminum. 
The  reason  being,  the  aluminum  shell  would  have  to  be  four  times  as 
thick  to  achieve  a  stress  field  of  equal  percentages  of  yield  stress. 
Hence,  titanium  will  always  give  a  longer  lifetime. 

While  the  sphere  furnishes  a  vastly  superior  creep  response 
than  the  cylinder,  it  is  more  expensive  to  construct  and  machine.   Thus, 
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the  engineer  must  weigh  the  merits  of  the  costs  of  materials,  con- 
struction and  machining  of  both  an  appropriate  sphere  or  cylinder  for 
his  particular  situation   (pressure  level  and  time  constraints)  before 
he  develops  a  final  design. 

Regardless,  the  shape  of  the  amplification  curves  for  both 
the  shell  shapes  suggest  a  creep  design  limit.   For  approximately  the 
first  seventy  percent  of  the  lifetime  imperfections  remain  small  and 
grow  in  a  linear  manner.   Thus,  an  arbitrary  design  limit  can  be  estab- 
lished which  states  that  as  long  as  the  imperfection  growth  becomes 
linear,  the  design  will  not  creep  buckle.   The  in-service  monitoring 
of  the  creep  deformation  process  is  easily  accomplished  with  a  three 
element  strain  gage  rossette  placed  at  the  point  of  greatest  imperfec- 
tion.  When  the  effective  stress  approaches  the  magnitude  of  the  yield 
stress,  the  active  use  of  the  structure  should  be  terminated. 
6. 5  Example  Calculations 
6.5.1   DSRV  ALVIN 

The  research  submersible,  ALVIN,  owned  by  the  Office  of 
Naval  Research  and  operated  by  Woods  Hole  Oceanographic  Institution, 
has  a  spherical  titanium  621/. 08  alloy  pressure  hull.   The  40.005  inch 
radius  hull  has  a  a/h  ratio  of  20  and  a  maximum  out-of-roundness  ratio 
(a*)  of  0.002.   The  maximum  operating  depth  is  12,000  feet  (5280  psi) 
and  a  collapse  depth  of  18,000  feet  (7920  psi). 

The  secondary  creep  law  for  Ti  6Al-2cb-lTa-0 . 8Mo  alloy  at 
room  temperature  has  been  empirically  measured  by  Chu  (1972)  (Appendix  C) 
While  the  measurements  were  done  in  tension,  it  is  believed  that  the 
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compressive  creep  behavior  is  similar  to  that  of  tension.   The  law  is 
expressed  as 

7       * 

o'    =  ke  (6.2a) 

where 

7 

k   =   8.3941934  x  104°  (-^-)   hr.        (6.2b) 

in2 
With  values  of  a/h ,  a*,  M  equal  to  20,  .002,  and  7  respectively, 

the  amplification  curve  of  Figure  34  is  manifested.   If  subjected  to 

the  maximum  operating  pressure,  buckling  occurs  at  a  time  equal  to 

9870  hours.   If  the  submersible  is  operated  at  only  half  its  operational 

depth  it  would  not  buckle  until  1.85  x  10°  hours  had  past. 

6.5.2   SOFAR  Floats 

The  V7oods  Hole  Oceanographic  Institution  has  employed  a  long 

cylinder  which  when  released  at  sea  sinks  to  a  depth  at  which  it  becomes 

neutrally  buoyant.   The  current  measuring  floats,  first  designed  by 

Swallow  (1963) ,  will  follow  a  water  mass  of  a  constant  density  due  to 

the  neutrally  compressible,  buoyant  properties.   Placed  at  a  depth  of 

approximately  5000  feet,  (the  axis  of  the  SOFAR  Channel,  a  horizontal 

channel,  in  which  the  convex  shape  of  the  sound  velocity-depth  profile 

allows  sound  to  be  transmitted  with  minimal  attenuation  over  thousands 

of  miles) ,  the  floats  emit  an  acoustic  tone  at  a  given  frequency. 

The  float  is  then  spatially  fixed  using  acoustic  triangulation  from 

land  based  listening  stations.   The  floats  have  been  successful  in 

tracking  eddies  in  the  WIIOI  MODE  program.   Unfortunately,  after  two  to 

three  months,  these  floats  have  been  found  to  sink  an  additional  2  30  to 

560  feet,  with  accelerating  sinking  velocity,  (an  apparent  creep  buckling 
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phenomenon).   Mavor  (1974)  computed  a  linear  strain  rate  of  0.9  uin/in/hr 
from  the  change  of  displacement  and  an  assumed  constant  in  situ  seawater 
density  gradient. 

-4rr  =  1.371  x  10~6/ft  (6.3) 

pAZ 

where 

p  =  density 

Z  =  depth 
The  floats  are  standard  ALCOA  extruded  6061-T6  aluminum 
cylindrical  stock  with  a  mean  radius  of  5.625  inches,  a  a/h  ratio  of 
7.5,  an  approximate  out-of- roundness  ratio  of  0.005  and  a  length  of 
17  feet.   They  operate  at  a  mean  depth  of  5000  feet  (2200  psi) .   The 
6061-T6  aluminum,  follows  a  cubic  creep  law  as  experimentally  found  by 
the  author  (Appendix  C) 

a3  =  kc  (6.4a) 


where 


k  =  4  x  1019  (psi) 3hr.  (6.4b) 


Figure  35  displays  the  predicted  creep  response  with  a/h, 
a*,  M  equal  to  7.5,  0.005,  and  3,  respectively. 

Since  the  problem  is  concerned  with  only  the  first  three 
months  of  the  lifetime,  the  perturbation  transverse  displacement  can  be 
considered  linear.   From  Figure  35 

w'  =  1.13795  yin/hr.  (6.5a) 

From  equation  (4.1) ,  the  perturbation  hoop  strain  rate  is 

e1'   =  ^™1  =  5.3947  x  10~8  in/in/hr.        (6.5b) 
°     a2 
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From  equation  (4.7b) ,  the  dominant  hoop  strain  rate  is 

•    _ 
—    w 

e_=  —    .  (6.5c) 

6    a 

Computing  w  from  equation  (4.8)  yields  a  dominant  hoop 
strain  rate  of 

e.  =  9.45  x  10~7  in/in/hr.  (6.5d) 

u 

Hence,  the  total  predicted  linear  strain  rate  is 

efl  =  eQ  +  e'   =  .998  yin/in/hr,  (6.5e) 

0     W       0 

ten  percent  greater  than  the  estimate  of  Mavor. 

This  calculation  shows  that,  in  fact,  creep  deformation  is 
most  probably  the  cause  of  the  unexpected  sinking  of  the  floats. 
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7.   SUMMARY  AND  CONCLUSIONS 

The  theoretical  investigation  which  has  been  undertaken  into 
the  creep  instability  of  the  infinitely  long  cylindrical  and  complete 
spherical  thin  shells  has  satisfied  the  original  thesis  objective  of 
developing  a  relatively  simple  analytical  method  to  predict  the  creep 
deformation  process  of  these  shells.   It  was  found  that  in  both  shells 
the  critical  buckling  mode  is  the  second  elliptical  mode.   The  roles  of 
shell  radius  to  thickness  ratio  and  initial  out-of-roundness  have  been 
explored.   The  relative  superiority  of  the  spherical  shell  over  the 
cylindrical  shell  has  been  discussed. 

The  investigation  demonstrated  that  the  single  most  important 
factor  in  creep  design  is  designing  the  structure  in  a  manner  which 
maintains  the  stress  field  at  the  lowest  possible  level,  thus  extending 
the  structure's  lifetime.   This  one  factor  endorses  the  use  of  high 
strength  titanium  and  the  spherical  shell  in  creep  design. 

In  conclusion,  while  creep  buckling  can  be  a  problem  in  long 
term  ocean  structures,  the  nature  of  the  time  history  of  imperfection 
growth,  almost  negligible  until  the  last  ten  percent  of  the  lifetime, 
allows  the  ocean  engineer  to  develop  a  safe  design  resistant  to  creep 
buckling  with  little  difficulty. 


8.   RECOMMENDATIONS 

Two  essential  advances  in  creep  design  are  required  in  the 
future.   First,  the  establishment  of  accurate  creep  constants  at 
ambient  temperatures  for  common  ocean  engineering  materials.   Second, 
a  more  advanced  theoretical  investigation  which  encompasses  the 
effects  of  boundary  conditions  of  such  shells  as  finite  length  cylinders, 
and  partial  spherical  shells  (end  caps) . 
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Figure  27       Non-dimensional  critical  time  with  various  initiation 
points  of  time  segments  of  constant  effective  stress 
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Figure  28       Critical  time  for  various  initial  elastic  loadings , 

expressed  as  a  percentage  of  yield  stress,  for  spheres 
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Figure  33       Ratio  of  critical  time  of  the  complete  sphere  to 

critical  time  of  the  cylinder  for  various  a/h  ratios 
(m  -  7). 

Figure  34       Creep  response  of  the  titanium  spherical  hull  of 
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Figure  35       Creep  response  of  the  aluminum  cylindrical  SOFAR 
float. 
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APPENDIX  A  -  COMPUTER  UTILIZATION 


(i)   Computer  Program  A  -  Perturbation  contour  map  of  the  complete 

sphere. 

The    following  FORTRAN  program  was   developed   to   compute 

w'/B^      f°r   a  discrete   number   of  points,    (5     by   5      grid)    in   a  given 
2m 

quadrant.      One   needs   only   to   input   the   initial   values   of   <j>    and   8 
of   a  given  quadrant   to   receive   the   x-y  position   and  value   of  w'/B^ 
of   the  before   mentioned  grid. 


x  -   sintj)cos0  (1) 

y   =    costj)  (2) 

2 
w  t— >       m 


^     P      (cos<}>)  cosm6  (3) 


B^  *-*         2 

2m        m=0 

PERTURBATION    CONTOUR   MAP    OF    THE    COMPLETE    SPHERE 
DIMENSION   T(19)  ,    P(19)  ,    W(19)  ,    E(19)  ,    F(19)  ,    P2(3) 
DO   51   M=l,    4 
READ(8,200)     U,V 
200    FORMAT    (2F10.6) 
DO   52    N=l,    19 
A=N-1 

T  (N)=.0174»5.0«A+U 
DO   53    1=1,    19 
B=I-1 

P(I)=.0174»5.0«B+V 
S=SIN(P(I)) 
C=C0S(P'(I)) 
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P2(l)  =  .5»(3.0oO»2.0-1.0) 

P2(2)=3.0»S»C 

P2(3)=3.0«S««2.0 

W(I)=0.0 

DO  54  J=l,  3 

R=J-1 

D=COS(T(N)»R) 

WZ=D»P2(J) 
54  W(I)=W(I)+WZ 

E(I)=S»COS(T(N)  ) 

F(I)=C 

WRITE(5,202)  E(I),  F(I),  W(I) 
202  FORMAT ( ' ■ ,  3F15.5) 
53  CONTINUE 
52  CONTINUE 
51  CONTINUE 

END 


Notation  of  Prog 

ram 

A: 

T 

e 

P 

+ 

P2(l) 

P  (COS0) 

P2(2) 

P   (cos<J>) 

P2(3) 

2  , 

P   (cost})) 

W 

transvers 

E 

X 

transverse  displacement/amplification  function 
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F  y 

U  initial  6  value 

V  initial  <f)  value 

(ii)   Computer  Program  B  -  Growth  of  the  amplification  function  of  the 

infinitely  long  cylinder. 

This  program  was  devised  so  that  one  inputs  a/h ,  a*,  M  and 

\\)  /ty     and  equation  (4.19)  is  computed.   The  output  consists  of  non- 
2   1 

dimensional  time  (t  ) ,  normalized  amplification  function  (A  /C  )  and 

c  2    2 

a  second  non-dimensional   time  which   is    a  better   indicator   of   real 

'tp/ 
time 


INFINITELY    LONG    CYLINDER 

MS  =  3 

JS=15 

DO    1    M=l,    MS 

READ(8,100)     B,D,E,XM 

100  FORMAT (4F10. 6) 
T=0.0 

Y=1.0 

WRITE  (5,101)     B,D,E,XM 

101  FORMAT (' 1" ,5X, ' A/H= ' ,F10.6, 'A  = ' ,F10.6,10X, 'PHI-' ,F10.6,10X, • M 
1='  ,F10.6,//,15X, 'TIME'  ,20X,  'AMP    FN '  , 10X ,  '  REAL   TIME',//) 

DO    2    J=1,JS 

IF    (J.GT.l)     GO   TO    7 
F=1.0 

GO   TO   8 

7    Y=E»X 

F=Y««(XM-1.0) 
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8  S=EXP(T»F/3.0) 

A=E»(1.0+6.0«E«B»D*EXP(T«F/3.0) ) 

T=3.0/F»ALOG( ( A-l. 0) / (B»D»F»6 . 0)  ) 

Tl=T/(9.0  B»»(XM+2.0)»(.866«(l.+6.  B  D)  )•  »(XM-1.)) 

WRITE(5,103)  T,A,T1 
103  FORMAT (' ' ,7X ,F10 . 3 ,17X,F10 . 3 ,15X ,E20 . 5) 
2  X=Y 

1  CONTINUE 
END 

The  number  of  input  statements  is  controlled  by  the  variable 
MS,  while  the  number  of  cycles  for  a  given  set  of  inputs  is  controlled 
by  the  variable  JS . 


Notation  of  Program 

B: 

T 

T 
C 

A 

A  /C 

2   2 

Tl 

tP  m 
z 

k 

Y 

* 

XM 

M 

B 

a/h 

D 

a* 

E 

2   1 

(iii)  Computer 

Program 

C 

-  Gro 

cation  function  with  time, 
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This   program  is   structured  exactly   as   the   cylindrical  use. 

Equation    (5.1    )     is    corrrouted    from   the    inputs    R/H,    R*   M   and  Tp  /ty    . 

1      2 

The   output   received  is   non-dimensional  time    (t    ) ,    amplification    function 

(B^    /D^    )    and   a  second  non-dimensional   time  which   is   a  better  indicator 
2m     2m  ,       mx 

of   real  time! — - —  1  .      Again,    the  number  of   input   statements   and  number 
of   cycles    for   a  given  set  of   input   is    controlled  by   the   variables   MS, 
JS    respectively. 


COMPLETE    SPHERE 

MS=8 

JS=50 

DO    1    M=1,MS 

READ(8,100)     B,D,E,XM 

100  FORMAT (4F10. 6) 
T=0.0 

Y=1.0 

WRITE  (5, 101)    B,D,E,XM 

101  FORMAT (*1'  ,5X, 'A/H='  ,Fl0.6, 'A=   '  ,F10.6,10X, ,PHI='  ,F10 .6 ,10X , *M 

1=' ,F10.6,//,15X, 'TIME' ,20X, 'AMP    FN ' , 10X, ' REAL   TIME',//) 

DO   2    J=1,JS 

IF    (J.GT.     1)     GO   TO    7 

F=1.0 

GO   TO    8 

7  y=E»x 

F=Y»»(XM~1.0) 

8  S-EXP(T»F/3.0) 
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Rl=4.5»B«D»F*S 

R2=112.5»(B»D»F»S)»  «2.0 

G=E  •  «2  .  0 • (0  .  5+R1+R2 )  -0 .  5 

A=((-R1/S)  +  ((R1/S)»   e2.0+4.0«G«»R2/S<»«2.0)«  »0  .  5)  /  (2  .  0«R2/S«  »2  .  0) 

T=3.0/F«ALOG(A) 

Tl=T»2./(9.«B«»(XM+2.)«  (l./2.«  •.  5»(  .5-4.  5»B»D+112.5»  (B«D)»  »2  .  )  ) 

0.5»(XM-1.) 

WRITE  (5  ,103)    T,A,T1 
103    FORMAT C ' ,7X,FlO . 3 , 17X,F10 . 3 , 15X,E20 . 5) 
2    X=Y 
1    CONTINUE 

END 


Notation  of  P 

rog 

ram 

C: 

T 

T 

s 

Y 

* 

A 

2m  2m 

Tl 

tP  m 

8 

k 

XM 

M 

B 

R/H 

D 

R* 

E 

2   1 

-106- 


APPENDIX  B  -  EQUIVALENT  NOTATION  FROM  PAST  INVESTIGATIONS 

(i)   Hoff,  N.J.  ,  Jahsman,  W.E.  and  Nachbar ,  W.  "A  Study  of  Creep 
Collapse  of  a  Long  Circular  Cylindrical  Shell  Under  Uniform 
External  Pressure",  J.  of  Aero  Space  Science,  Vol.  20,  No.  10, 
1959. 

Notation  of  Paper  Equivalent  Notation 


aR 


7-  =  (^)(!f)V)-2 

3   AA    p 


1  9    9 

t    =  ^  log(l  +  T/X  2) 
cr   2         2   o 


2a*  (£, 


I(Pam  2a 2 
kv  h'  (  h1 

o„i  + 1  ft  4 

8    a   a* 


2t 
h/1 


_k_ 
,  n 


K 


n  M 

where   p  =  —  radius  of  gyration 
R  =  average  radius 
a  =  shape  factor 
X  =  weighted  shape  factor 
(ii)   Berman,  I.,  Chern,  J.M.  and  Gupta,  G.D.  "A  Parametric  Study 

of  Elastic-Plastic  Creep  Buckling  of  a  Thin  Cylindrical  Shell", 
J.  of  Pressure  Vessel  Technology,  August  1974. 
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Notation  of  Paper  Equivalent  Notation 


+ 
,.    w  -w- 

0 


0    2a 


q  P 

(iii)   Samuels on,  L.A.  "Creep  Buckling  of  a  Cylindrical  Shell  Under 

Non-Uniform  External  Loads",  Int.  J.  Solids  Structures,  Vol.  6, 
1970. 

Notation  of  Paper  Equivalent  Notation 

cf>  6 

B  K 

R  a 

The  remainder  of  the  notation  is  consistent. 
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APPENDIX  C  -  CREEP  CONSTANTS  OF  COMMON 
OCEAN  ENGINEERING  MATERIALS 

(i)   Titanium  621/. 08  (Ti-6AL-lTa-. 8  Mo) 

(a   =  120  ksi) 

y 

Chu  (1972)  found  that  the  creep  data  of  Ti  621/. 08  at  room 
temperature  can  be  represented  as 

,       o        .13. 890  0. 183 
E  =  {I£^)  t  (1) 

where 

e  =  creep  strain,  % 
a  =  stress ,  ksi 
t  =  time ,  hours 
If  equation  (1)  is  plotted  at  a  given  stress  level,  it  is 
discovered  that  after  an  initial  period  of  primary  creep,  the  non- 
linear viscous  behavior  of  secondary  creep  becomes  apparent 

as  example:   for  o  =  103.55  ksi,  e   =  .1071  yin/in/hr. 
By  calculating  the  constant  strain  rates  for  a  number  of 
stress  levels  the  following  power  law  emerges 

o7  =  8.3  x  10k°l  (2) 

where  a  =  stress,  psi 

e  =  creep  strain  rate,  yin/in/hr 
(ii)   Titanium  6-4  (6AL-4V) 

(a  =  135  ksi) 

y 

Chu  has  found  the  emperical  equation  for  titanium  6-4  at 
room  temperature,  but  has  not  as  yet  published  it. 
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(iii)   Aluminum  6061-T6  (wrought  alloy  -  magnesium  and  silicon  alloy 

solution  treated  and  artifically  aged) 
(c   =  40  ksi) 

y 

The  author  conducted  a  number  of  tests  in  which  cylindrical 
columns  of  Aluminum  6061-T6  were  loaded  with  a  uniform  load  in  compres- 
sion by  a  lever  system.   Deformation  was  measured  using  two  uniaxial 
strain  gages  mounted  on  opposite  sides  of  the  column  in  a  series  con- 
figuration.  The  series  configuration  allowed  bending  effects  to  be 
internally  cancelled  out  of  the  readings  by  the  gages. 

The  Aluminum  6061-T6  was  found  to  abhor  to  a  cubic  law  at 
room  temperature 

a3  =  4  x  1019e  (3) 

where 

a  =  stress,  psi 

e  =  creep  strain  rate,  yin/in/hr. 
All  tests  were  conducted  at  high  load  levels,  greater  than 
65  percent  of  yield,  at  room  temperature.   The  strain  rates  computed 
are 

a  (psi)  e(uin/in/hr) 

34670  0.8504 

33750  0.8420 

31450  0.7768 

27615  0.0569 

(iv)   Stainless  Steel  304  (Austentic) 

Berman  (1974)  expresses  the  secondary  creep  behavior  of 
Stainless  Steel  304  at  1050°F  as 
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where 


a5-6789  _  8%035  x  lo28e 

a  =  stress,  psi 

e   =   creep   strain,    yin/in/hr. 
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